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ABSTRACT 

We investigate the conditions under which the regular satelhtes of Jupiter and Saturn formed. 
The final stage of giant planet accretion is thought to occur slowly over a relatively long, 10 Myr, 
timescale. Gas accretion during this stage, through a completely or partially opened gap in the 
solar nebula, occurs slowly allowing for the condensation of ices, and incomplete differentiation, 
seen in the regular satellites of the giant planets. Furthermore, the dichotomy seen in the Jovian 
and Saturnian systems may be explained as this infall wanes or is completely shutoff as a result of 
gap opening or global depletion of gas in the solar nebula. We present one-dimensional simulations 
of circumplanetary disks that couple the viscous transport of material with the loss of mass at 
the disk outer edge by ultraviolet photoevaporation as well as the infall of material from the 
solar nebula. We find that the circumplanetary disks of these protoplanets are truncated, as a 
result of photoevaporation, at a range of values with the mean corresponding to ~ 0.16 Hill radii. 
These truncation radii are broadly consistent with the current locations of the regular satellite 
systems of Jupiter and Saturn. We also find that photoevaporation can successfully act as a 
clearing mechanism for circumplanetary nebulae on the potentially short timescales, 10^ — 10^ 
yr, over which mass accretion from the solar nebula wanes as a result of gap opening. Such a 
rapid clearing of the circum- Jovian disk may be required to explain the survival of the Galilean 
satellites. 



Subject headings: planets and satellites: general 

1. Introduction 

In recent years, the impact of photoevapora- 
tion on the evolution of p rotoplanetary disks has 
received much a t tention (Ijohnstone et al.l 1998t 



solar system: formation 



the evolution of protoplanetary nebulae ( Clarkd 



20071: iMitchell fc Steward 120101 ). It is tempting 



ive 

20091 ). Recent observations confirm that there is 
sufficient FUV and X-ray flux from young, solar- 
type stars to drive photoevap orative mass lo ss in 
the surrounding nebulae (Jnglebv et al. I lionh . Al- 
though photoevaporation has been investigated in 
detail in the context of planet formation, it has 
yet to be applied to circumplanetary disks and the 
formation of regular satellites. 

Recent models have looked at the combined ef- 
fects of viscous evolution and photoevaporation on 



Ballv et al. 1998: Adams et al. 



20041 : 



Gorti et al 



to apply a direct analogy between protoplanetary 
and protosatellite nebulae, but it would be in- 
correct to view a protosatellite nebula as simply 
a scaled down version of a protoplanetary neb- 
ula. The slow infall of material from the solar 
ne bula onto circumplanetary d is ks, as proposed 



bv ICanuD fc WardI (|200a l2009l) : IWard fc Canupl 
POIOI ). causes them to behave very differently 
from protoplanetary analogues. In particular, the 
combination of mass infall and photoevaporative 
mass loss can allow for quasi steady-state solu- 
tions. 

The giant planets have similar ratios of the total 
mass contained within the satellite systems to the 
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mass of the host planets (Msateiiites ~ 10 Mpia„ct) 
In general, the mass of a satellite increases with in- 
creasing distance from the planet, reaches a maxi- 
mum, and then decreases with distance. A striking 
difference is seen in the mass distributions of the 
satellite systems of Jupiter and Saturn. Jupiter 
has four large moons of relatively the same size 
whereas Saturn has one large moon that contains 
almost all of the mass in the entire satellite sys- 
tem. The similarities in these systems indicate 
that they formed from similar processes, yet these 
processes have acted in such a way to produce 
very different outcomes. 

The low inclination and prograde orbits of 
the regular satellites of giant planets indicate 
that they formed in situ from a circumplane- 
tary disk. The similarity of the satellite sys- 
tems of giant planets, with Jupiter in particu- 
lar, to the solar system has contributed to the 
development of formation th eories that employ a 
"min imum mass subnebula" (jLunine &: Stevenson 
19821 ) , analogous to the "m inimum mass solar neb- 
ula" ( Weidenschilling|[l977 ) . Summing the mass of 
heavy elements contained within the satellites and 
augmenting it to solar composition determines the 
mass of this subnebula. The "MMSN" approach 
leads to a minimum mass of ^ 0.02 Mj. This how- 
ever leads to a variety of problems when applied 
to the Jovian and Saturnian systems. 

For one, the resulting nebula is too warm to 
condense ices, yet, with the exception of lo, con- 
densed volatiles are known to be a major compo- 
nent of the Galilean satellites. In such a dense 
nebula, the effect of type I migration would be 
strong enough to migrate Callisto into Jupiter be- 
fore it had suffic ient time to grow to its current 
size (|Warcil 1 19981 ). Even if type I migration did 
not operate, Callisto would accrete on such a short 
timescale it would melt and become differentiated, 
which is not support ed by Galileo observations 
(jAnderson et al.l 119971 ). Recent gravity measure- 
ments by the Cassini spacecraft indicate that Ti- 
tan is also partially differentiated which would im- 
ply a long accretion timesc ale for the regular satel- 
lites of Saturn as well (|Iess et al. 201(A Barr et al. 
l201fl[ ). 

Beginning with the assumption that not all 
of the 0.02 Mj, estimated in the "MMSN" ap- 
proac h, be present at once in the protosatellite 
disk, ICanup fc WardI (|2002l ) have developed a 



model referred to as the "gas-starved" disk model. 
Their gas-starved disk has a slow infall of mate- 
rial, from which the satellites are slowly built. 
This slow infall results in a steady-state disk with 
a much lower instantaneous disk mass than in the 
aforementioned "MMSN" approach. In the "gas- 
starved" disk model satellites preferentially grow 
in the inflow region of the disk where gas and 
solids are delivered from the solar nebula. Once 
the satellites become massive enough, they mi- 
grate toward the planet, sweeping up solids along 
the way. The growth and loss by migration of suc- 
cessive generations of satellites allows for a quasi 
steady-state in which the total mass contained in 
the satellite system is maintained at ~ 1 0~'^ Mj 



(ICanuD fc WardI 120061: ISasaki et al.ll2010l) . When 
the disk mass is depleted as a result of the global 
depletion of the solar nebula, the radial migration 
of satellites is halted, only the final generation of 
satellites is retained. 

Furthermore, the observed ice-to-rock fraction 
of the Jovian satellites can be reproduced from the 
radial temperature gradient of the circumplane- 
tary disk in Canup & Ward's (2002) model. In the 
context of this model, Callisto forms on a longer 
timescale and thereby also avoids differentiation. 
The rapid inward migration problem is also reme- 
died in the "gas-starved" disk model. This model 
has also been applied to the Saturnian system to 
successfu lly explain the in complete differentiation 



of Titan (jBarr et al.ll20ia ) 



A slow infall rate, and resultant low disk mass, 
is also consistent with current theories of the for- 



mation of Saturn's i ce-rich ring system (jCanup 



2010l ). ICanupl (|2010l ) asserts that Saturn's ring 



system was formed when a differentiated. Titan- 
sized satellite migrated to within the Roche limit 
and its icy mantle was tidally disrupted. The 
migration of such a satellite requires that some 
amount of remnant gas still be present in the 
circum-Saturnian nebula. An upper limit to the 
amount of gas present is set by the requirement 
that the ring system not be removed before the 
circumplanetary nebula completely disperses. 

Recently, IWard fc Can"^ (|2ni(l[ ) have devel- 
oped a comprehensive model that follows the for- 
mation and evolution of a giant planet, and sub- 
sequent circumplanetary nebula, from cloud col- 
lapse through the contraction phase. Their model 
is composed of three elements: 1) an inflow model 
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describing the properties of the in-flowing material 
from the circumstellar nebula, 2) a quasi steady- 
state disk model, and 3) a planet growth and con- 
traction model. It is their quasi steady-state vis- 
cous disk model that most concerns us. 

The quasi steady-state viscous disk model 
requires an in-plane flux as well as a mass 
loss mechanism at th e outer disk boundary, 
(IWard &: Canud l20ld) . Despite the comprehen - 



Ward fc Canut 



gp renen - 



sive nature of their model, 
have yet to identify an appropriate mass loss mech- 
anism and in their own words, "the outer edge of 
the disk, r^, is not well defined other than it be 
much further out than the centrifugal radius". 
Given a moderate external FUV flux from either 
the central star or nearby high-mass stars, photoe- 
vaporation could provide a natural mechanism for 
both mass loss and truncation at the outer edge 
of circumplanetary disks. 



CanuD fc WaTdl (|2002h assume the same cir- 



cumplanetary disk evolution for both Jupiter and 
Saturn and that the dichotomy present in their 
satellite systems is explained by the stochastic 
timing of formatio n/migration and the depletion 
of the solar nebula ( Canup &: Wardl[2006l) . In con- 
trast, Sasaki et al. ( 201Cl[ ) surmise that the differ- 
ence in the satellite systems of these two planets is 
the result of very different disk evolution scenar- 
ios. Numerical simulations show that Jupiter is 
large enough to open a complete gap in the solar 
nebula cutting off infall onto t he Jovian subnebul a 
on a 10^ — 10"* year timescale ( Sasaki et al. 2010f) . 



In contrast, Saturn is too small to have opened 
a complete gap in the solar nebula and infall onto 
the Saturnian subnebula would have halted on the 
~ 10® year timescale for the dissipation of the solar 
nebula as a whole. Although the assumption that 
Jupiter opened a gap and Saturn did not depends 
on the assumed viscosity and scale height of the 
protoplanetary disk, the critical mass for opening 
a gap is larger in the outer regions of the disk. 
It is a reasonable assumption using standard disk 
parameters, and further supported by the current 
masses of Ju piter and Saturn. For these reasons, 
Sasaki et al.l (|2010, ) assume in their models that 
Jupiter opened a gap in the solar nebula and Sat- 
urn did not. 

Sasaki et al. I (|2010l) have recently published re- 
sults from a suite of simulations in which the 
growth and dynamical evolution of proto-satellite 



embryos was modeled. They seek to explain why 
the Jovian regular satellite system consists of four 
nearly equal mass satellites whereas the Saturnian 
system contains one large satellite. They propose 
that the Jovian satellite system may have been 
"frozen" in place when Jupiter grew sufficiently 
large that a gap was opened in the solar nebula. 
In the Saturnian system, where only an incomplete 
gap may have formed, the slower shutoff timescale 
for material infalling from the solar nebula would 
have allowed its satellites to continue to dynami- 
cally evolve. The typical end result is significant 
depletion of solids in the inner disk and the re- 
tention of a single large satellite in the outer disk 
that is similar to Titan. This outcome may also 
be a natural result of the inside-out clearing of the 
solar nebula due to photoevaporation by the Sun 
(see discussion in Section |3] below). 

Sasaki et"all (|2010[ ) are able to produce four or 



five similarly sized satellites in the Jovian system 
in 80% of their runs. Whereas, in the Saturnian 
system only one large satellite remains in 70% of 
their runs. Their models however rely on a rapid 
dispersal mechanism for circumplanetary gas once 
the gap in the solar nebula has been opened. Pho- 
toevaporation could provide just such a mecha- 
nism as well as help to de termine a natur a l oute r 
disk boundary As with Ward k Canud ( 2010l ). 
the models of ISasaki et aL ~ (j201o[ l relv on an ad 
hoc outer disk boundary. 

Recent simulations have investig ated the tidal 
trunc ation of circumplanetary disks ([Martin fc Lubow 
20101) . These tidal truncation simulations produce 
disks that are truncated at a radius that occurs at 
^ 0.4 rjj, where ru is the planet Hill radius. 



I Mp 
3Mc 



(1) 



where a and Mp are the planet's semi-major axis 
and mass and Mq is the mass of the Sun. This 
outer disk radius is, however, too large to explain 
the compact configuration of the regular satellite 
systems of Jupiter and Saturn which extend to less 
than 0.06 rn. 

The primary objective of the present paper is 
to show how a nominal FUV fiux can photoevap- 
orate the out e r port ions of circumplanetary disks. 
Adams et al.l (|2004l ) demonstrate how photoevap- 
oration creates a subsonic outflow of gas in the 
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disk that is well inside the gravitational radius, 
rg, where the thermal velocity of the hot disk at- 
mosphere equals the planet's e scape velocity. We 
apply the lAdams et al. photoevaporation 
model to circumplanetary disks and find that the 
disks are truncated well inside the gravitational 
radius, rg. It is important to note that photoevap- 
oration only needs to remove gas from the planet's 
Hill sphere in order to truncate the circumplane- 
tary disk; it does not need to remove the gas from 
the solar system. 

The difference in solar flux, due to the dif- 
ference in semi-major axes, may further account 
for the different evolutionary histories of Jupiter 
and Saturn. Saturn's greater distance from 
the Sun would cause the incident solar flux to 
be a factor of (aj/as)^ less at Saturn than at 
Jupiter, where oj and as are the semi-major 
axes of Jupiter and Saturn respectively. Even 
in the compact configuration proposed in the 
Nice model, the difference in incident flux would 
be (5.45 AU/8.I8 AU)^ = 0.44, implying that 
the incident flux would have been more than 
twice as strong at Jupiter's location than Sa turn's 



(|Tsiganis et al.ll2005l: iMorbidelh et al.ll2005l) . The 



greater amount of incident flux at Jupiter would 
have caused there to be a greater amount of pho- 
toevaporative mass loss in the Jovian system. Fur- 
thermore, the increased rate of photoevaporation 
would have caused the Jovian subnebula to be 
more drastically truncated than the Saturnian 
subnebula given the same solar luminosity. 

We describe our disk model as well as our mod- 
els for infall from the solar nebula and photoevap- 
oration in Section[2l The results of our simulations 
are presented in Section [31 with the results from 
our steady-state disk models presented in Section 
13.11 and the results from our time-dependent, de- 
caying disk model in Section [321 In Section [3] we 
present a summary and discussion of our results. 

2. Disk Model 



Our disk model is an extension of an earlier pho- 
toevaporative viscous disk model that was previ- 
ously applied to the evolu tion of the solar nebula 
(iMitcheh fc Stewartll2010[) . Besides its application 
to circumplanetary nebulae, the current model dif- 
fers from the previous model in that it includes in- 
fall from the solar nebula. The details of the infall 



will be discussed in Section 12.11 It is this infall 
that allows for steady-state solutions. 

Our model uses the common a-viscosity pre- 
scription with a viscosity that is proportional to 
r. 

V = J^o(r/i?o) (2) 

where and i?o are scalings for viscosity and ra- 
dius. This assump tion has been used in the past 



by m any authors (jClarkd l2007t iHartmann et al 
19981 ). The linear dependence of viscosity on ra- 



dius in our model implies that the temperature 
profile is proportional to r~^l'^ . The midplane disk 
temperature, Tdisk7 is used in order to determine 
the viscosity constant vq. 



IkT, 



disk 



(3) 



where k is the Boltzmann constant and (/z) is the 
mean molecular weight. We have scaled the mid- 
plane temperature profiles in our models such that 
in the Jovian system the temperature has been 
set to « 250 K at 10 This is consistent with 
the slow-inflow, low-opacity circu m- Jovian accre- 
tion disk model investigated by (|Canup fc Ward 
20021 ). For the Saturnian system, the radial tem- 
perature profiles in our models have been set to 
« 100 K at 20 i?s, the approximate location of 
Titan. 

Most models of circumplanetary disks employ 
a value of a that is generally s maller than that 
used in circumstellar disk models (jCanup fc Ward! 
20091) . The simulations presented in this work con- 
sider a values in the range of 10~^ — 10~^, with 
our reference models having a value of a = 10"'^. 

Instead of using r and E; the radius and mass 
surface density of the gas, it is useful to describe 
the system in terms the variables h and g; the 
specific angular momentum and torque. This is a 
good choice of variables to make when the viscos- 
ity is proportional to r as it allows us to trans- 
form the viscous disk equation into a simple, lin- 
ear differential equation with a constant coefficient 
(jHartmann . 19981. 



,1/2 



dr 



(4) 



(5) 
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where G is the gravitational constant and f2 is the 
Keplerian orbital velocity. The relationship be- 
tween g and S in Eqn. ([5]) implies that the mass 
surface density is inversely proportional to the vis- 
cosity. Thus, E (X . 

By substituting in the functional form of viscos- 
ity, the continuity equation for radial mass trans- 
port can be written as 



dg i^oGMp d^g 
dt " 4 Rn dh^' 



(6) 



4 Ro 

where Mp is the mass of the planet. 

As we are interested in the addition of material 
to the disk, a source term must be added onto the 
right-hand side of Eqn. 



dg 3 i^oGMp d'^g Stti/q 



dt 4 i?o GMpRa 



h't{h,t), (7) 



where t^{h,t) is the infall rate in [g cni^^s^^] as 
a function of specific angular momentum (radius) 
and time. 

2.1. Infall 

The partial differentiation of Callisto implies a 
lon g formation tim cscalc and was one factor that 
led ICanup fc Ward (,2002) to develop their "gas- 
starved" disk model. They assert that the 0.02 
Jupiter masses required for formation the Galilean 
satellites, using a minimum mass solar nebula ap- 
proach, need not be present all at once. This ma- 
terial may have slowly flowed through the system 
over an extended period of time. In this scenario, 
the agglomeration of solids into satellites is analo- 
gous to the build-up of mineral deposits in plumb- 
ing over time. They introduce a timescale for the 
addition of mass to the system. 



TG 



Mp 

Mn 



(8) 



In all simulations presented here, we use a gas ac- 
cretion timescale of tg = 5x10^ yr. This is equiva- 
lent to the mass accr etion rate of 2 x 10~^M p yr~^ 
which was found by ICanup fc WardI (l2002l) to be 
consistent with the conditions needed to form 
the Galilean satellites. This mass accretion rate 
is also consistent w ith gap accretion simulations 
(jBrvden et al.lll999l ). 

The location of the outer boundary is not only 
governed by photoevaporation, but by viscous 



spreading as well. It is therefore important to con- 
sider the location and amount of infalling material 
from the solar nebula. If the majority of the in- 
falling material falls onto the outer regions of the 
disk, because of its angular momentum content, 
it will affect the behavior of the outer boundary. 
It is therefore important to parameterize the lo- 
cation and amount of infalling material, and not 
.ius t treat it as a constan t thro ughout the disk as 
do lCanup fc WardI (|2002l l2006l) . 

Machida et al. I (l2008l) ran high resolution three- 



dimensional simulations of the accretion of angu- 
lar momentum onto a protosatellite system. Us- 
ing their results, we can estimate the centrifugal 
radius, rc, the location at which the bulk of the 
infalling gas lands on the circumplanetary nebula. 
We estimate that, for a Jupiter mass protoplanet, 
the infalling material has an average specific angu- 
lar momentum consistent with a Keplerian orbit at 
Tc « 25 We use the angular momentum con- 
tained within the infalling material to constrain 
the radii where infall occurs, so that the added 
disk mass has the same angular momentum than 
the local circular orbit where it is added. Were this 
not the case, and the infalling material contained 
less angular momentum as the local circular orbit 
where it was added, it would not be rotationally 
supported and would rapidly fall inward and be 
redistributed at smaller radii corresponding to its 
angular momentum content. A very recent, high- 
resolution simulation of mass accreted from the 
solar nebula onto circumplanetary disks indicates 
that infalling mass does, if fact, intercept the disk 
at radii (~ 50i?j) greater than that which corre- 
sponds to its angular momentum content and is re- 
distributed inward (Tanigawa et al. 2011,). While 
these results are too new to have been included in 
this investigation, they will certainly be taken into 
consideration in future investigations. 

We have adopted the following functional form, 
which peaks at 25 rp, for the infalling mass onto 
our protosatellite nebulae. 



t{r,t) = - 







"(^) - 28' 




cxp 
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g cm ^s ^, 



(9) 

where rp is the planet's current equatorial ra- 
dius. The prescription for the infalling material 
has roughly 90% of the infalling mass accreting 
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onto the circumplanetary disk at radii between 
16 i?p < r < 28 Rp and peaks at 25 Rp. This 
prescription is c onsist ent with the simulations of 



Machida et al. 



(|2008l) which show that the in- 
fahing material has some angular momentum dis- 
tribution around the peak which corresponds to a 
centrifugal radius of, Vc — 25 Rp. Our prescrip- 
tion for infall is also c onsistent with the previous 
work of ICanup fc War d (2002, 2006) in which the 
infall is limited to an inner region of the disk. In 
our model there is no accretion of material onto 
the circumplanetary disks from the solar nebula 
external to r = 28 i?p. 

Once the gas giant planets grow large enough, 
they begin to open gaps in the solar nebula as a 
result of resonant interactions with the disk. The 
gap-clearing timescale can be estimated by assum- 
ing that it would occur on the viscous timescale to 
spread ac ross the scale heig ht in the protoplane- 



tary disk (jSasaki et al.ll201 



r 



Tgap ^ ^ ^ (10-3 _ 10-4) X (1 - 10) Myr (10) 



The actual timescale for a growing planet to open 
a gap is likely longer than this estimate, but it 
must be substantially shorter than the planet's ac- 
cretion timescale in order to limit the final mass 
of the planet. A reasonable estimate for the gap 
opening timescale would be to assume a median 



value of Tc 



gap 



2.5 X lO'^ yr. 



We assume that the infall rate decays exponen- 
tially over some timescale. Toff. 



S(t) = t exp 



t 

Toff 



(11) 



This assumed form for the infall decay rate is also 
made by oth er authors in similar investigations 
(ICanuDll2O10l) . We further assume that the infall 
decay time is of the same order as the gap opening 
timescale, Tos — Tgap — 2.5 x 10^ yr. 

Although it seems likely that gas continues to 
accrete through the gap in a circumstellar disk 
opened by a growing protoplanet, it is uncer- 
tain if this gas is able to accrete either onto the 
protoplanet or even onto a circumplanetary disk. 
We have chosen the exponential decay of infalling 
material shown in Eqn. [TT] because it is simple 
and easy to unde rstand. Also , it all ows a direct 
comparison with ISasaki et al. ()2010l) . which use 



the same prescription for the exponential decay 
and short infall decay timescales. Furthermore, 
whether it be via gap opening or the global deple- 
tion of the solar nebula, it is certain that the infall 
from the solar nebula should halt at some point in 
time. The only uncertainty is the timescale over 
which the infall wanes. 

2.2. Photoevaporation 

Recent observations of 56 weak-lined and clas- 
sical T Tauri stars was used to determine the FUV 
emission of young 1 — 10 Myr stars and whether 
the observed flux is consistent with that required 
by models of ph otoevaporative disp ersal of cir- 



cumstellar disks ( Inglebv et al. 201ll) . These au- 



thors concluded that radiation fields sufficiently 
strong for the removal of gas are present during 
the disk dispersal phase. Although not considered 
in this study, x-ray flux has beei i shown to enhanc e 
FUV-driven photoev aporation ( Gorti et al. 2009() . 
Inglebv et al.l ( 201l[ ) also investigated the x-ray 
emission of 1 — 10 Myr stars and found that the 
x-ray flux remains high, and constant, throughout 
the duration of the dispersal phase. 

Clearly, the young Sun would have exposed the 
circumplanetary nebulae of Jupiter and Saturn to 
FUV radiation. Although it is unclear at this time 
how much FUV from the Sun is able to reach 
these circumplanetary nebulae, one must remem- 
ber that the Sun most likely formed in a cluster of 
10'^ — 10^ stars which could have contributed to the 
FUV field in which these nebulae wer e embedded 
(jFatuzzo fc Adamsll2008l : lAdamsll2O10l ) . Whatever 
its source, the FUV radiation would have heated 
the periphery of the circumplanetary disks. Gas 
heated to sufficient temperatures would then have 
become unbound from these disks. The gravita- 
tional radius, rg, is defined as the radius at which 
the sound speed of the heated gas equals the es- 
cape speed from the system. 



kT 



(12) 



where k is the Boltzmann constant and T is the 
temperature of the super-heated atmosphere, or 
what we will refer to as the envelope temperature, 
Tonv Gas beyond the gravitational radius will es- 
cape from the system. The gravitational radius 
is the canonical radius beyond which gas heated 
to a temperature Tcnv will escape from the disk. 
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In actuality, gas can escape from the disk at radii 
substantially smaller than rg. 

The heating of the disk's surface by FUV radia- 
tion and the resultant outflow are complicated pro- 
cesses, but it is useful to employ a simplified model 
with an isothermal, heated atmosphere. Consider 
a disk irradiated and heated by external FUV ra- 
diation from either the early Sun or the fellow 
members of its birth cluster. Depending on the 
strength of the FUV flux, the heated gas will reach 
temperatures in the range 100 K < T < 3000 K 
(Adams et al.. 2004). As the gas heats, it expands 
generating a neutral outflow. The expanding out- 
flow begins subsonically but becomes supersonic 
by the time it reaches the gravitational radius. 
This outflow is generally isotropic, but the ma- 
jority of mass loss is dominated by mass loss from 
the outer edge of the disk. The isotropic, neu- 
tral outflow serves to shield the disk from EUV 
radiation that would io nize the disk an d heat 
it to » 10,000 K See Ijohnstone etall ([1998, ): 



Adams et al.l ()2004l ) for more thorough discus- 



sions. 



Adams et al. ( 20041 ) studied the previously un- 



explored subcritical regime, where the outer radius 
of the disk, r^, is smaller than the gravitational 
radius. Using results from their photodissociation 
region code, they developed analytical models for 
the photoevaporative mass loss rates for cases in 
which rd, the location of the outer edge of the 
disk, is both inside and outside the gravitational 
radius. These models are characterized by a sin- 
gle temperature, Tenv, which determines a unique 
sound speed, Og, for the isothermal atmosphere. 

M.CoiVc(.>a...(^)exp(-^; 

Td < /3rg (13a) 

M = 47rJ"(^)crpuyasrd 

rd > /3rg (13b) 

The first equation is for subcritical disks, and the 
second equation is for supercritical disks, where 
as is the sound speed in the heated atmosphere 
and Co is a constant of order unity. Nq is the 
critical surface density of the flow and ctfuv is 
the cross section for dust grains interacting with 
FUV radiation. The dust optical depth is given 
by Tpuv = CTpuv ■ For an optical depth of 

order unity, cTpuy ~ -^h, where A'h is evaluated 



at the critical density A^c ^ 10 cm~ . The factor 
Co is a constant of order unity used by Adams et 
al. (2004) to match their numerical and analytical 
solutions. It is used in our model to match the 
mass loss rates for sub- and supercritical disks at 
a radius of ^d/rg — 0.25 and is adopted as our 
value of /3. This is necessary because the mass loss 
rates are sensitive functions of strength of the FUV 
radiation field as well as the assumed matching 
point. 

We have matched the subcritical solution onto 
the supercritical solution, because the supercrit- 
ical solution is better understood and well con- 
strained. The factor J^, in the second equation, 
is the fraction of the solid angle subtended by the 
flow and is ^ 1 because the flow from the disk sur- 
face and edge merge at a radius between and 
2 • rd: creating a nearly s pherically symmetric out- 
flow (jAdams et al.ll2004l ). 

The differential equation governing the lo- 
cation of the outer boundary was derived us- 
i ng mass conservation a t the outer boundary 
(jMitchell fc Stewartllioiflh . 



Mb 



= M 



Af photo- 
cvaporatic 



(14) 



The mass flux due to viscous processes is simply 
—dg/dh and the mass flux due to photoevapora- 
tion is taken to be Eqn. (|13aP or Eqn. (jl3bp de- 
pending on the location of rd. The mass flux due 
to the motion of the boundary is iirr^^.T'didrd/ dt). 
Here, and throughout the paper, the subscript "d" 
indicates that these quantities are to be evaluated 
at the outer disk edge. The equations governing 
the location of the outer boundary were derived 
by substituting these expressions into Eqn. (|14p . 
Depending on the location of the outer boundary, 
we will either be in the subcritical regime. 



dr| 
dt 



2ttNc{ij,) dhd 
or the supercritical regime 



exp 



dt 



dg 



2nNc{fi) dhd 



2tT QsTd, 



(15) 



(16) 



At each time step, depending on whether we 
are in the sub- or supercritical regime, either Eqn. 
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((15)1 or Eqn. p6|) must first be solved to determine 
the new boundary location. Then the viscous evo- 
lution, Eqn. ([7]) is solved. All of the simulations 
presented here were performed on a Mac G5 Pow- 
erPC. They were evolved on a grid of 201 points 
evenly spaced in non-dimensional specific angular 
momentum. 

Our initi al conditions are that of the si milarity 
solutions of iLvnden-Bell fc Pringlel(ll974l) . 



Mq 

2TTRir 



exp 



r 



(17) 



where Ri is the initial disk scaling radius and A/q 
as the initial disk mass. 

As with most previous disk models, we employ 
a zero torque inner boundary condition. This al- 
lows for accretion from the disk onto the planet. 
Our outer boundary condition is set by the torque 
exerted on the disk by the out-flowing material 
that is being photoevaporated at the disk's edge. 
The torque exerted on the outer edge of the disk 
can be expressed as 



5d = 3V2^^Nc{^I)lyo 



Rq 



(18) 



It must be noted that these two boundary con- 
ditions are used to solve the equation governing 
the temporal evolution of the mass surface density. 
We have, in effect, a second outer boundary condi- 
tion that governs the temporal evolution of the lo- 
cation of the outer boundary. This condition is set 
by a differential equation that was derived using 
mass conservation at the outer boundary (Eqn.'s 

(usi) & (uni)). 

2.3. Variable Space Grid Method 

A well-posed problem in the material sciences, 
called the Stefan Problem, deals with propagat- 
ing phase changes typically considered in t he con- 
text of melting/freezing and heat ablation (jOzigik 
19801) . We have developed a one-dimensional 



model of an astrophysical disk that includes vis- 
cous diffusion and photoevaporation at the outer 
boundary using numerical te chniques developed 
to so lve the Stefan Problem (jMitchell &: Stewart 
l2010[ ). By adapting the Stefan problem to as- 
trophysical disks, our one-dimensional numerical 
model self-consistently tracks the location of the 



outer boundary. This is a novel approach to mod- 
eling astrophysical disks. 



To solve the Stefan problem, iKutluay et al 



(|1997[ ) adopt a numerical method with a variable 
space grid (VSG) first proposed bv lMurrav fc LandisI 
( 1959t ). The VSG method employs a fixed number 
of grid points with a variable grid size at each time 
step. This method involves solving two coupled 
differential equations at each time step, one for 
the location of the outer boundary and one for 
the diffusive evolution of the disk. Once the loca- 
tion of the outer boundary is found the abscissa 
is rescaled and the diffusive evolut i on ca lculated. 
Please refer to iMitchell fc Steward (|2010[ ) for fur- 
ther details. 

3. Results 

3.1. Steady-State Disks 

We present a number of simulations of steady- 
state circumplanetary disks about Jupiter and 
Saturn. These simulations produce steady-state 
solutions by coupling viscous evolution, photoe- 
vaporation and infall from the solar nebula. A 
goal of this work is to investigate the range of 
possible truncation radii and disk masses based 
on the amount of FUV fiux and the strength of 
the viscosity parameter, a. 

These simulations all agree well with t he low 
surface densities found b ylCanup &: Ward (120021) . 
As in the simulations of ( Machida et al. 2008|) . we 



see surface density enhancements at roughly 25 Vp, 
as a result of infall from the solar nebula. Such en- 
hancements are not seen in recent, high-resolution 
simulations in which the infall occurs over a wide 
range of radii (jTanigawa et al.l 120111) . These en- 



hancements likely arise from our choice of infall 
shape (Eqn. ([9|)) and the fact that our viscosity 
depends only on radius and not on the local condi- 
tions in the circumplanetary disk. If such density 
enhancements are real, they may have implications 
for satellite growth. We will further explore the 
significance of these density enhancements in Sec- 
tion H 

The UV flux can produce a wide range of en- 
velope temperatures (100 K — 300 K) depend 



ing on the magnitude of the flux ([Adams et al 



20041 : iMitcheh fc Stewa"rtll2010l ). As stated in Sec- 
tion l2.21 recent observations confirm that the FUV 
emission of young, solar-type stars is sufficient 
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for photoevaporation ( Inglebv et al.|[201ll) . Fig- 
ure [1] shows the radial mass surface density from 
three simulations in which the envelope temper- 
ature has been varied. The solid line is of our 
reference model with Tenv = 600 K, whereas the 
dotted and dashed curves are for simulations with 
Tenv = 100 K and Tenv = 3000 K respectively. 
These simulation were all run with a = 10"'^. One 
significant feature of these simulations is the en- 
hancement in mass at r 25 Rj mentioned above. 
The truncation radii of these simulations ranges 
from 73 Rj in our high-temperature simulation 
to 324 Rj in our low temperature simulation. In 
terms of Hill's radius, these disk radii range from 
0.098 — 0.44 rn with our fiducial model's outer ra- 
dius at 0.17 th. 




100 200 
radius [Jupiter radii] 

Fig. 1. — Steady-state radial mass surface den- 
sity of Jupiter's circumplanetary disk. The three 
curves are for three different values of the isother- 
mal, heated envelope. The solid line is of our 
reference model with Tenv = 600 K, whereas the 
dotted and dashed curves are for simulations with 
Tenv = 100 K and Tenv = 3000 K respectively. The 
current location of Callisto has been marked with 
a dot-dashed line. 



The higher envelope temperature causes more 
erosion at the outer boundary and therefore causes 
the steady-state disk to have a much smaller trun- 
cation radius. Despite having a much smaller 
truncation radius, the three runs presented in Fig- 
ure[T]all have disk masses that lie in a narrow range 
from 3.3 x 10"^ - 1.2 x lO"'' Mj. 

We have also run a suite of simulations in which 
the viscosity parameter, a, has been varied. We 



have investigated a in the range of 10"'* — 10~^. 
Simulations of steady-state circum- Jovian disks 
are presented in Figure [2 The solid line is of 
our reference model with a = 10^'^, whereas the 
dotted and dashed curves are for simulations with 
a = 10~* and a = 10~^ respectively. These simu- 
lation were all run with Tenv = 600 K. A striking 
difference between the varied viscosity runs and 
those for a varied envelope temperature, presented 
in Figure [U is that the outer disk radius is inde- 
pendent of the strength of the viscosity. It depends 
only on the mass loss rate at the outer boundary 
which, in our models, is controlled solely by the 
envelope temperature. 




50 100 
radius [Jupiter Radii] 

Fig. 2. — Steady-state radial mass surface den- 
sity of Jupiter's circumplanetary disk. The three 
curves are for three different values of the viscos- 
ity parameter, a. The solid line is of our reference 
model with a = 10~^, whereas the dotted and 
dashed curves are for simulations with a = 10"'* 
and a = 10~^ respectively. The current location of 
Callisto has been marked with a dot-dashed line. 



Another difference between the temperature 
and viscosity runs is in the masses of the steady- 
state disks. The disk masses in the runs presented 
in Figure [2] each differ by an order of magnitude 
and lie in the range of 5.5 x 10"^ - 5.4 x 10~* Mj. 
This suggests that the steady-state disk mass is 
dominated by the strength of the viscosity and 
not by the amount of photoevaporative mass loss. 

We have conducted a similar suite of simula- 
tions for the Saturnian system as for the Jovian 
system. Simulations in which we have varied the 
envelope temperature are presented in Figure [3] 
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and those in which we have varied the strength 
of the viscosity are presented in Figure ID The 
same range of parameter space has been explored 
in these runs as was done for the Jovian sys- 
tem. These data are presented with the same line 
scheme as those presented in Figures [T] and [2] 



Saturnian Circumplanetary Disk 

variiUion in viscosity pammclcr cc 




50 100 
radius [Saturnian radii] 



150 



Fig. 3. — Steady-state radial mass surface den- 
sity of Saturn's circumplanetary disk. The three 
curves are for three different values of the isother- 
mal, heated envelope. The solid line is of our 
reference model with Tenv = 600 K, whereas the 
dotted and dashed curves are for simulations with 



100 K and To, 



3000 K respectively. The 



current location of Titan has been marked with a 
dot-dashed line. 



Due to its location farther out in the solar 
nebula, the conditions during formation were 
probably somewhat different at Saturn than at 
Jupiter. The inside-out clearing of th e solar neb- 
ula by photoevaporation described by iGorti et al.l 
pOOgf ) would have caused the solar FUV flux to 



vary significantly in space and time. It is how- 
ever almost certain that both Jupiter and Sat- 
urn experienced some migra t ion d uring the for- 



mation process (jWardl Il997l |1998[) . The direc- 
tion and rate of migration are highly uncertain 
and we cannot s ay for sure where either of these 
systems formed (^ Paardck oooer fc Mellemal 2006t 
Crida fc Morbidelli.,2007, : , Walsh et al.ll2nill ). Fur- 



thermore, we feel that our range of parameter 
space is appropriate for both the Jovian and Sat- 
urnian systems and exploring the same range in 
the two systems is best for comparison. 

Once again, one can see in the simulations of 




20 40 
radius [Saturnian radii] 



60 



Fig. 4. — Steady-state radial mass surface den- 
sity of Saturn's circumplanetary disk. The three 
curves are for three different values of the viscos- 
ity parameter, a. The solid line is of our reference 
model with a — 10~^, whereas the dotted and 
dashed curves are for simulations with a = 10~^ 
and a = 10~^ respectively. The current location 
of Titan has been marked with a dot-dashed line. 



the Saturnian system that the total disk mass is 
dependent on the strength of the viscosity whereas 
the radial extent of the modeled disks is dominated 
by the envelope temperature. The disk masses 
for each disk presented in Figure |4] also vary by 
an order of magnitude, just like those presented 
in Figure [2] and lie in the range of 2.3 x 10^^ — 
1.9 X 10"* Ms. The outer disk radii of the disk 
models presented in Figure [3] span a range from 
38 — 137 Rs, which in terms of the planet's Hill 
radius is 0.035 — 0.13 th. Despite the wide range 
of outer disk radii, the simulations presented in 
Figure [3] have disk masses in a narrow range from 
1.0 X 10^5 - 5.1 X 10"^ Afs. 

The time-dependent, viscous evolution of our 
models allow us to accurately track the transfer of 
mass throughout our disks. The transfer of mass is 
important for the formation of satellites. The rate 
of mass transfer through the disk, M can be calcu- 
lated at any radius by simply taking the derivative 
if the torque with respect to the specific angular 
momentum. 

d£ 
dh 



M = 



(19) 



Disk models are often characterized by the 
slope of the radial mass surface density using a 
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power law of the form 



E(r) cx 



(20) 



The slope of the radial surface mass density pro- 
files of the models presented in Figures [2] and |4] 
are roughly power laws with q = —1.3. This slope 
is ste eper than the slopes u s ed in the dis k mod - 
els of ICanuD fc WardI (l2002l ): ISasaki et alJ (l2010l) . 
who both assume q = —3/4. However, it is in 
accordance with our assumed radial temperature 
dependence Of r^^/^, which predicts q ~ —1. The 
fact that q is slightly steeper than q = —1 in 
our models can be attributed to the truncation 
of the outer boundary by photoevaporation. It 
must be kept in mind that all of the slopes men- 
tioned here are ad hoc and have been assumed by 
the various authors. It is unclear how steep these 
disks would actually be given a more comprehen- 
sive model with realistic viscosity that depends on 
the local conditions in the disk. 

Figure [S] shows the mass transfer rate as a 
function of radius in a steady-state circum- Jovian 
disk. This analysis was conducted on our refer- 
ence model, with an envelope temperature of 600 
K and a viscosity parameter a — 1Q~^. Nearly 
equal mass is transferred inward and outward in 
the disk, with slightly more mass accreted onto 
the planet than is lost through the outer edge due 
to photoevaporation. Of the 2 x 10^^ il/j yr~^ 
of material being accreted from the solar nebula, 
1.25 X 10~^ Afj yr~^ is accreted onto Jupiter and 
0.75 X 10^^ il/j yr^^ is lost through photoevap- 
oration at the disk's outer edge. Assuming a so- 
lar dust-to-gas mass ratio of 0.014, this is suffi- 
cient to provide enough mass in solids through the 
outer regions of the disk to build Callisto over the 
10^ yr required for it to remain undifferentiated 
(ICanuD fc Wardll2002[ ). 



Aft, 



1.9 Afcallisto 



(21) 



where / is the assumed dust-to-gas mass ratio of 
0.014 and race is the accretion timescale of 10^ yr. 

3.2. Decaying Infall 

A secondary goal of this work is to validate pho- 
toevaporation as a potential mechanism for the 
rapid dispersal of circumplanetary nebulae as in- 
fall from the solar nebula wanes because of gap 
opening. We have performed one such simulation 




40 60 80 100 120 
radius [Jupiter radii] 

Fig. 5. — Mass transfer rate in a steady-state 
circum- Jovian disk as a function of radius. This 
analysis was conducted on our reference model, 
with an envelope temperature of 600 K and a vis- 
cosity parameter a — 10^'^. Nearly equal mass is 
transferred inward and outward in the disk, with 
slightly more mass accreted onto the Jupiter than 
is lost through the outer edge due to photoevapo- 
ration. The mass accretion rate in the outer region 
of the disk is characterized by a constant mass flux 
rate. The region of constant flux begins exterior 
to 28 Rj. Interior to 28 Rj, the mass flux rate 
is decreases and at ^ 24 i?j it becomes negative. 
The region inward of 28 i?j corresponds to the re- 
gion over which infall from the solar nebula occurs. 
Inward of ~ 24 Rj the mass flux rate is negative 
and is being accreted onto the planet. 



that investigates the decrease in infall rate onto 
the Jovian subnebula over a 2.5 x 10^ yr timescale 
(see Section 12. 1|) . Figure |6] shows the temporal 
evolution of Jupiter's circumplanetary disk as the 
infall is abated. 

Similarly to lSasaki et al. I ()2010l) . we have begun 
this simulatio n in the context of the "gas-starved" 
disk model of lCanup fc WardI (|2002l ). This is done 
because gap opening would likely occur in the fl- 
nal stage of giant planet accretion. Therefore, the 
simulation presented in Figure [S] was begun with 
the steady-state solution from our Jovian reference 
model with tq = 5 x 10^ yr. This value difi^ers 



slightly from that used by ISasaki et al. ( 2010l ) of 



2 x 10^ yr, but as they have shown the distribution 
and composition of final satellites is insensitive to 
the exact value of tq. The radial surface density 
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50 100 
radius [Jupiter radii] 

Fig. 6. — Temporal evolution of Jupiter's circum- 
planetary disk as the infall from the solar neb- 
ula exponentially decays. The radial mass surface 
density is shown at 500yr increments, with the 
solid bold lines indicating t — Toff, 2 -Toff and 3 -Toff. 



is shown at 500 yr increments, with the solid bold 
lines indicating t — Toft,2 • Toff and 3 • Toff. The 
total mass of the circumplanetary disk decreased 
nearly two orders of magnitude from 5.5x10"^ Mj 
to 7.6 X 10~^ Af J over the course of this simulation. 

This simulation was done using a nominal FUV 
flux, which corresponds to an envelope tempera- 
ture of 600 K. Envelop e temperatures can range 
from 100 K to 3000 K (| Adams et al.l 12004 ). This 
shows that, even for a moderate FUV flux, photo- 
evaporation is able to clear the Jovian subnebula 
on the very short timescale over which the infall 
wanes due to gap opening. 



4. Summary and Discussion 

We have modeled and analyzed circumplane- 
tary, protosatellite disks with the combined in- 
fluence of viscous forces and photoevaporation. 
Our models also include mass infall from the solar 
nebula allowing for steady-state solutions. These 
models were developed in the con text of the "gas- 
starv ed" disk models put forth by (jCanup fc Ward 
20021 ). With these models, we present a new mech- 
anism for the truncation of circumplanetary disks. 
We model the late stages of giant planet growth, 
when the accretion rate from the solar nebula is 
limited. The limited accretion rate may be a re- 
sult of gap opening, the global depletion of the 
solar nebula or a combination of both. 



The regular satellites of Jupiter and Saturn 
extend to 0.06 th- This value is much smaller 
than the extent of the circumstellar disk pre- 
dicted considering the angular momentum content 
of the accreting gas. By estimating the angu- 
lar momentum of accreting gas as it travels from 
the solar nebula, through the Lagr a nge p oints, 
and onto Jupiter, iQuillen fc Trillin esti- 
mated that the truncation radius is ~ rH/3. A 
similar outer disk radius was found in numerical 
simulations which consider only t he angular mo- 
ment um content of accreting gas (lAvliffe fc Bate 
2009 ). The rece n t, hig h- resolution simulation 



of iTanigawa et al.l (|2011l ) indicates that accret- 
ing gas intercepts the Jovian circumstellar disk at 
~ ^h/IS and that much of the gas has an angu- 
lar momentum content that corresponds to even 
smaller radii than this. 

Recently, the radial extent of the regular satel- 
lites of Jupiter and Saturn has been explained 
by the rapid accumulation of solids at the loca- 
tion of the infalling material, Vc or by the trun- 
cation of the p r otosa te llite disk by solar tide s 
(ICanup fc WardI 120021 : iMartin fc Lubo^ HoIS). 
We invoke photoevaporation to naturally explain 
the location of regular satellites. Our models pro- 
duce a large range of outer disk radii depending 
on our choice of envelope temperature (FUV flux) . 
These outer disk edges range in radii from 0.035 rn 
to 0.44 rn, with a mean value of 0.16 rn. Even 
though our models show a wide range of outer 
disk edge radii based on envelope temperature, 
the masses of these disk all lie in a narrow range 
from 3.3 x 10~^ - 1.2 x 10"" Mj in the Jovian 
models and 9.7 x 10-^-5.1 x 10"^ Ms in the Sat- 
urnian system. This narrow range of disk masses 
is to be expected because these simulations were 
all run with the same alpha-viscosity parameter 
of 10-3 

The disks produced in our Saturnian models ex- 
tend to a much smaller extent of their Hill spheres 
because Saturn is less than a third as massive as 
Jupiter and nearly twice as far from the Sun. In 
relation to their masses, Saturn's Hill sphere is 
much larger than Jupiter's. Although the size of 
a planet's Hill radius plays no role in our simula- 
tions, similarly sized disks, relative to their host 
planet's mass, would extend to different fractions 
of their respective Hill spheres. 

Our models show that photoevaporation can 
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truncate circumplanetary disks to radii that are 
in concordance with the locations of the regular 
sateUites of the giant planets in our solar sys- 
tem. Our reference model produces outer disk 
radii truncated at 0.057 th and 0.17 rn for Sat- 
urn and Jupiter respectively. These small outer 
disk radii may provide an obvious explanation for 
the locations of the regular satellites of Jupiter 
and Saturn. For the outer disk radii to set the 
scale of the present regular satellite systems, the 
solids from which the satellites formed must have 
been small enough to be coupled to the gas disk. 
This could occur if small solids are continually be- 
ing supplied from the solar nebula or through col- 
lisional processes within the proto-satellite disk. 
However, if the solids grow on a timescale that 
is much shorter than the timescale for orbital de- 
cay due to gas drag, then the solids will decouple 
from the gas and the radial extent of the resulting 
satellite system would be set by the centrifugal ra- 
dius, Tc. The location and transport of solids in an 
actively supplied, photoevaporating circumplane- 
tary disk is an important issue that merits further 
investigation. 

The constant outward mass flux in the outer re- 
gions of the circum- Jovian disk, presented in Fig- 
ure [5l indicate that while very low surface densi- 
ties exist in the outer regions of these disks there is 
still significant amounts of gas being transported 
to this region (see Figure [5|). Assuming that solids 
are carried along with the gas as it is transported 
outward, there is sufficient mass present for satel- 
lite formation. Given a solar abundance of solids, 
we calculated that there is nearly twice as much 
mass in solids transported outward over a 10^ yr 
time period than is needed to form Callisto. There 
are reasons to believe that the solids-to-gas mass 
ratio would be higher than solar and therefore we 
take this value as an underestimate. We plan to 
investigate the actual distribution and transport 
of solids in the near future with a more compre- 
hensive model that is currently in development. 

Another aspect of our simulations that may 
have consequences on satellite formation is the 
density enhancement seen in our simulations 
at ~ 25 Tp, the location of peak mass in- 
fall. A similar den sity enhancement was seen in 
Machida et al.l ( 20081 , but was not present in more 
recent simulations wh ich include radiative transfer 
(lAvliffe fc Batell2009[ ). It is unclear at this time if 



these density enhancements would be present in 
a more realistic model in which the viscosity was 
determined locally. This is an interesting question 
and one which we plan to investigate in the near 
future. If it is real, these enhancements would 
have a significant impact on satellite formation. 
Density enhancements such as these are accom- 
panied by pressure maxima. It has been shown 
that migrating solids can be trapped in such pres- 
sure maxima and ra pidly grow into satellitesimals 
(jKretke et al.ll2009l) . 

In an effort to test what effect the location at 
which infalling material intersects the disk has on 
steady-state disk morphology, we have performed 
one test simulation in which the peak of the in- 
falling material occurred at 35 Tp rather than at 
25 Tp. In the test simulation, the location of the 
disk outer edge was shifted farther out by ~ 3% 
and the total disk mass increased by ^ 5%. The 
changes are a result of a greater fraction of the 
infalling mass being transported outward rather 
than inward. This test indicates that it is im- 
portant to identify the exact location at which 
the infalling mass accretes onto the circumplane- 
tary disk. However, the current 3-D hydrodynam- 
ical simulations used to model infall from the so- 
lar nebula onto circumplanetary disks have insuf- 
ficient resolution to identify the location exactly. 
We are satisfied that our treatment is sufficient for 
this study, yet we plan on including more precise 
results as they become available. 

While the strength of the viscosity plays no role 
in the location of a disks outer boundary, it does 
play a significant role in the total mass contained 
in a given disk. In a steady state, the mass accre- 
tion rate, M oc ^S, is constant. A constant mass 
accretion rate implies that the mass surface den- 
sity must be proportional to the inverse of the vis- 
cosity. In our models, this means that the surface 
density is inversely proportional to the viscosity 
parameter, a. One would naively assume that a 
larger surface density would result in larger satel- 
lites, but in actuality the opposite is true due to 
the increased rate of migration. A more massive, 
lower viscos ity disk results i n a les s massive satel- 
lite system. Canup &: WardI (2002)find that satel- 
lites will only survive against type I migration for 
values of a > 10""^, assuming a solar gas-to-solid 
ratio for the infalling material. 

By varying the viscosity parameter, a, we are 
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able to produce a wide variety of total disk masses 
in our models. These masses range from 2.3 x 
10~^ — 1.9 X 10"'* Ms in the Saturnian system 
to 5.5 X 10"'' - 5.4 X 10""* A/j in the Jovian sys- 
tem. These radii are all many of orders of magni- 
tude smaller than the ^ 0.02 A/j inferred for the 
"MMSN" approach. Such low mass surface den- 
sities are a result of the 5 x 10^ yr timescale over 
which the infall occurs. Again, these conditions 
are necessary to produce the ice-rich compositions 
of the regular satellites of Jupiter and Saturn as 
well as the incomplete differentiation of Callisto 
and Titan. 

The location of the disk outer boundary in our 
models depends only on the amount of FUV flux, 
as parameterized by our envelope temperatures. 
These radii are independent of the viscosity. This 
may not however be the case in a more realistic 
case where the viscosity is dependent on the local 
properties in the disk and is not simply propor- 
tional to the radius. We are currently developing 
a more sophisticated model in which the viscosity 
will be calculated locally. We seek to explore the 
morphology of steady-state disks with more real- 
istic viscosities in the near future. 

As discussed earlier, Jupiter, and to some ex- 
tent, Saturn, is expected to open a gap in the 
solar nebula as a result of resonant interactions. 
Whether the gap opened be complete or partial, it 
would seriously inhibit the flow of gas onto circum- 
planetary, satellite forming disks. The simulation 
presented in Figure [6] shows that if such a restric- 
tion in infall is accompanied by a rapid disper- 
sal mechanism, such as photoevaporation, it will 
cause the rapid removal of the circumplanetary 
disk as well. The rapid removal of gas would dy- 
namically "freeze" the Galilean satellites in place. 
The rapid dispersal may also be aided by the in- 
creased ionization fraction, and increased viscos- 
ity, that may occur as a result of lower surface 
densities. We hope to investigate this in the future 
with a model that includes the local evaluation of 
viscosity. 

The partial, or incomplete, gap opened by Sat- 
urn woul d potentially have very different conse- 

Again, photoevapo- 



the incomplete differentia tion of Titan and the for- 
mation of Saturn's rings (iBarr et al.ll2010l : ICanupl 
2010l ). We plan a more thorough investigation of 



quences (jSasaki et al. 



J nave 
20101) . 



ration would limit the size of the Saturnian sub- 
nebula based on the waning infall rate. The very 
low surface densities that would be present as a re- 
sult of the waning inflow would be consistent with 



the decay of infall in the near future. 

We feel that the gap-opening scenario described 
in Sasaki et al. 1 dkoio) can be further refined by 



including the effect of the inside-out clearing of 
the solar nebula that occurs in many simulations 
which include irradiation from a central source. 
In particular, we would hke to draw the re ader's 
attention to Figure 4 in lGorti et"al] (|2009[ ). The 
snap-shots of the mass surface density in Figure 
4 show the photoevaporation front sweeping out- 
ward past Jupiter's location rather rapidly, but 
taking much longer to sweep past Saturn's loca- 
tion. Although it is difficult to determine from 
the snap-shots shown in Figure 4, the outward- 
traveling photoevaporation front is moving at ap- 
proximately 8.8 AU Myr"^ through the region 
where Jupiter is located. This is roughly twice 
as fast as it is moving through the region in 
which Saturn is located, which is approximately 
4.5 AU Myr~^. The outward-traveling front not 
only sweeps past Saturn at a slower rate than 
Jupiter, but nearly 1 Myr later as well. 

The inside-out clearing of the solar nebula may 
have exposed Jupiter's circumplanetary disk to 
much more solar fiux at an earlier time than Sat- 
urn. Also, at Saturn's location the surface den- 
sity seems to be decreasing globally on a simi- 
lar timescale over which the front moves outward. 
If so, it may mean that the shutoff of the infall 
from the solar nebula onto the Jovian subnebula, 
because of the local depletion of the solar neb- 
ula due to inside-out clearing, would have hap- 
pened sooner and on a shorter timescale than in 
the Saturnian system. I n our scenario, the slow 
infall rates required by ICanuT) fc WardI (|2002h : 
Ward fc Can^ (|201(j) would occur as a result of 



gap opening, or partial gap opening in the solar 
nebula. Slow infall from the solar nebula, onto 
the circumplanetary nebulae, would continue un- 
til the outward-traveling photoevaporation front 
passes their respective locations, globally clearing 
the solar nebula and ceasing any further infall. 
In this scenario, the final structure of the regu- 
lar satellite systems of the giant planets that we 
see today would have been determined by the rate 
and timing of the inside-out clearing of the solar 
nebula. 
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On a final note, the processes responsible for the 
transport and removal of angular momentum in 
astrophysical disks is highly uncertain. It is even 
more uncertain in the case of circumstellar and cir- 
cumplanetary disks where the ionization fraction 
of the gas may be low enough to prevent magne- 
torotational instabilities from occurring. What- 
ever the mechanism for angular momentum trans- 
port, it must have occurred in circumplanetary 
disks to enable the accretion of the planet. Gas 
liberated through photoevaporation carries angu- 
lar momentum away with it and may thereby pro- 
vide a natural explanation for its removal from the 
system. 

This work was supported by NASA Earth and 
Space Science Fellowship #10-PlanetlOF-0057 
and by the Cassini Project. G.R. Stewart was 
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